---
title: "R9.IAN_Points"
output: pdf_document
---

# Load Packages
```{r}
library(tidyverse)
library(stargazer)
```

# Load Data
```{r}
load("Data/ian_points.Rdata")
```

# Run Regressions (App. Table 3ci)
```{r}
m.1 <- lm(log(tax_income) ~ log(points_total), data = ian_points)
m.2 <- lm(log(tax_income) ~ log(points_total) + log(area_total), data = ian_points)
m.3 <- lm(log(tax_income) ~ log(points_total) + log(point_ratio), data = ian_points)

stargazer(m.1, m.2, m.3, type = "text", 
          covariate.labels = c("Log(Total Points)", "Log(Total Area)", "Log(Irrigated Points Ratio)"), 
          dep.var.labels = c("Log(Taxable Income)"))
```

